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ABSTRACT 


Understanding the concepts based on problem solving is not an easy 
methodology in teaching the impact of power systems stabilizer (PSS) on 
transient synchronous generator using MATLAB capability. Experiments 
conducted in simulating sessions play an important role in this teaching. This 
simulation can simulate power system stability behavior with reasonable 
accuracy in less time. This transient phenomenon of a power system utilizing 
synchronous generator and modelling by fully three-phase model with 
changes in stator flux linkages neglected is analyzed by employed single 
machine infinite bus taken to the power system. Whereas a power system 
stabilizer which consist of a wash-out circuit, two stages of compensation, a 
filter unit, and a limiter, is applied to control voltage and frequency of power 
systems in transient condition. Proportional-integral-derivative (PID) 
controller tuned by Ziegler-Nichols’s method is cascaded to conventional PSS 
in order to enhance the response time of system while providing a better result 


in damping for oscillation. This gives the clear idea about PSS and PID 
controller impacts on transient synchronous generator and its enhancement to 
the students of electrical engineering program, Institut Teknologi Nasional 
Yogyakarta. 
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1. INTRODUCTION 

Teaching the power systems stability concepts regarding transient synchronous generator becomes a 
challenge. This begins from the fact that the power systems must preserve voltage and frequency in the 
required level under any disturbances, mostly electromechanical low frequency oscillations (LFO) [1], [2]. 
Under transient conditions, the voltages or the currents are still significantly sinusoidal in shape, but the 
amplitude fluctuates with time and the signal may be shifted away from the zero line [3], [4]. These 
oscillations decrease power transfer capability and can trigger instability in the power system, which in the 
end can lead to failure of the power system [5]. The one of many methods to guarantee stability is to damp 
electromechanical oscillations by using generators’ excitation or automatic voltage regulator (AVR) [6]. 
High gain excitation system in supplement with AVR will damp out the oscillations of low frequency with 
small amplitude. Another method is augmenting supplementary control loop to the input of the AVR and 
excitation system of synchronous generator called by the PSS or power system stabilizer [7]. PSS were 
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introduced as a mean to boost damping through the modulation of synchronous generator’s excitation so as to 
extend the power transfer limit. Input signals to the PSS may be proportional to angular rotor speed, 
generator output frequency or active power output. Indeed, the calculations involved in this power systems 
stability are very tiresome and complex. 

Understanding the concepts based on problem solving is not an easy methodology in teaching 
transient stability of synchronous generator in power system. Experiment plays an important role in such 
power system stability teaching. It has to affiliate firmly with theory teaching. The entering of computer 
simulation to the teaching process has become an inevitable trend [8], [9]. This simulation which possible to 
simulate power system stability behavior with reasonable accuracy in less time, will become an interface 
between the theory teaching and experimental session. Many electrical engineering departments teach power 
system stability based on simulation, even for the laboratory studies, but some number of these simulation- 
based studies fail to give the comprehensive concept of the system. There are many non-linearities in power 
systems stability, especially in synchronous generator models. These non-linearities cannot be completely 
brought out by the simulation tools due to its difficulty in implementation and time consumed by the 
computers. 

The non-linearities in power systems stability must be changed into a linearized model. Although 
modern control methods have been used to minimize the arranged objective function of power system 
stabilizer, namely fuzzy logic PSS [10] and adaptive fuzzy rule-based PSS [11], the conventional lead-lag 
structure of power system stabilizer is still chosen because of the ease to tune online. The reason behind that 
might be the ease to tune online and the lack of stability guarantee regarding to some methods of adaptive or 
variable structure [12]-[15]. Maintaining the best damping performance when there is an extreme change in 
system operating condition, an auto-tuning stabilizer of proportional-integral-derivative (PID) structure can 
overcome this problem. Many studies related to such structures have been carried out [16]-[18]. But, the PID 
parameters are still selected and optimized by adding various search algorithms, such as bacteria foraging, 
particle swarm optimization, and fuzzy logic controller. Those PID’s structures still leave a problem, namely 
the controller's stability guarantee [19], [20]; On the contrary, the PID controller tuned by Ziegler-Nichols’s 
method preserves the stability margin [21]. 

The simulator of conventional lead-lag model of PSS cascaded with PID controller for controlling 
power systems stability on transient synchronous generator has been developed and implemented in 
Electrical Engineering Program Study, Faculty of Industrial Technology, Institut Teknologi Nasional 
Yogyakarta. This simulator preserves the stability condition of PID structure tuned by Ziegler-Nichols’s 
method. 

This paper explains the way in which simulator of the power system stabilizer for stabilizing 
transient condition of synchronous generator is developed. Procedure for making instability problem and how 
the problem is to be corrected is demonstrated. Section 2 explains research method. Result and analysis are 
given in section 3. The conclusion is drawn in section 4. 


2. RESEARCH METHOD 
2.1. Transient stability of synchronous generators in power systems 

In modern power system where a number of synchronous generators are operated in parallel, studies 
are usually achieved to confirm that the generators will operate appropriately in the event of possible faults or 
changes in system conditions. Studies concerned with transient stability of synchronous generator examine 
the ability of the generators to preserve synchronism from large oscillation formed by a transient disturbance. 
Because of a large oscillation, models of machines should describe the crucial nonlinearity in the frequency 
range between 1.0 and 5 Hz [22]. The dynamic behavior of synchronizing oscillation is influenced by system 
parameters and kind of control. 

In a large power system, it is not practical to depict each and every component in full detail. 
Depending on the range of response frequency counted to be important for the problem at hand, proper model 
of the required fidelity would be selected. The electromechanical oscillation frequency of synchronous 
generator in power system lies between 0.5 to 3 Hz. But, the transient time constant ranging from 0.5 to 
10 seconds is usually longer than the period of the electromechanical oscillation. 

To reduce the effort spent on establishing the model and on computation, the simplification that are 
obtainable by the separation in time scales of the different dynamic behaviors and that the harshness of a 
disturbance is usually attenuated as it propagates through the system; 1.e., the duration of the electrical 
transient of the network is relative short to the electromechanical dynamic of synchronous generator. So, a 
static representation of the network can be used where longer electromechanical oscillation are primarily of 
interest. The models chosen for generators on the network need not be the same. Some understanding of 
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basic dynamic behavior of synchronous generator in transient condition can be attained using a simple model 
of generator and network, known as single machine infinite bus. 

The single machine infinite bus with neglected stator transient or stator gdO flux linkages is shown 
in Figure 1. If voltage drop across external line is V, = (% + jx,)I, using transient model then the equations 
of d-axis and q-axis winding voltages without contribution of the kd and kg damper windings will be 
changed; the parameters of linearized synchronous generator model is shown in Table 1. The derived sets of 
linearized equations of the system are AV, AEG; AE;, and Ao. Usually, small resistance of stator winding, 7, is 
neglected. 





generator 


Figure 1. Single machine infinite bus system [23] 


Table 1. Transient generator model with changes in stator gdO flux linkages neglected [23] 
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Taking slight displacement considering values of the steady state operating designated by and supplementary 
subscript o, we will get 
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2.2. Power system stabilizer 
2.2.1. Concept of power system stabilizer 

Heffron and Phillips [24] developed single machine connected to an external grid using an 
electromechanical model of a synchronous generator with an excitation system. Then, De Mello and 
Concordia [25] used such model to describe the electrical oscillations and damping torque in an electrical 
system, shown in Figure 2. In this figure, the contribution of torque created by the stabilizer path is given by: 


ATEP 


PN Kz Et 
AWG 


K6 OE ref 


|Exc(jw,)| 


PSS,,(s)GEP(s) & P(s) GEP(s) = IGEP(jw,)| = K, (3) 


T 


The GEP(s) is a transfer function which defines the characteristics of the generator, the excitation 
system, and the power system, shown in Figure 3. The generators typically have a voltage regulator transient 
gain which results in satisfactory operation conditions. Slight phase lag of excitation system will relatively 
give high gain at small loads, but will intentionally reduce damping of rotor oscillations and may affect an 
instability. 
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Figure 3. Transfer function blocks of synchronous generator with exciter and PSS using slip speed [27] 


The value of 1/(K3Tp,) is assumed less than the cross-over frequency. The gain is also proportional 
to the parameter K, which denotes the result of a change in generator flux (E4) on torque; consequently, it is 
increased with generator loading and also increased as the power system becomes stronger. Quite the reverse, 
for situations where the voltage regulator loop crossover frequency is greater than the oscillation frequency 
of concern, the gain of GEP(s) is no longer proportional to regulator gain. It is inversely proportional to the 
parameter K,, which decrease as the ac system becomes stronger and hence causes the gain of GEP(s) to 
further increase as the power system strength increases. 


2.2.2. Structure of conventional power system stabilizer 

The conventional PSS consists of wash-out circuit, two stages of phase compensation, filter and 
limiter. A wash-out circuit is a high pass filter (HPF) for readjust action to remove steady-state offset, shown 
in Figure 4. The two stages of phase compensation which composed of two lead-lag compensators have 


center frequencies of compensation of 1/27 ./T,T, and 1/27 ./T3T,. The lead stage is used to compensate 
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for the phase lag introduced by the AVR and the field circuit of the generator. Meanwhile, a filter section is 
used to suppress frequency components in the input signal of PSS that could stimulate undesirable interaction 
and a limiter prevents the PSS’s output signal from driving the excitation into saturation. The PSS’s output 
signal is fed an adding input signal, V;,,,,), to the regulator of excitation system. 

The AVR is to maintain the terminal voltage of generators because of adjusting the exciter voltage 
of the generators. AVR unit typically includes an error amplifier with limiters. The degree of transient gain 
reduction can be done using a compensator that has Tę < Tg. The stabilizer feedback signal, Vp, and 
supplementary signal, Vij), are arrived at the input regulator from a PSS. The AVR’s output must be 
increased by the exciter before it has the necessary power and range to excite the field winding of a large 
generator. The exciter consists of the field winding, magnetic nonlinearity of the exciter’s main field path, 
and the armature. The voltage equation of the field winding and region of S, curve around the normal 
Operation condition can be approach by an exponent function; its current and voltage as shown in (4). Figure 
5 shows a block diagram of the transfer function between v, and vf. The feedback gain is a function of rp, 


but Tg will not be affected by changes in 77. 





; : r dv 
lfpu = Venu Tr Rao on = |k; + i Vxpu + TE PE (4) 


where 







Figure 4. Power system stabilizer with AVR and exciter [28], [29] 
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Figure 5. Magnetization curve and transfer function of exciter (in pu) [30] 


2.3. Auto-tuning stabilizer of PID structure 

Since conventional PSS consists of dual lead-lag compensator unit and its parameters are not tuning 
by original system parameter, so its control is less flexibility and the control results are far from ideal during 
drastic higher changing of its operation condition, such as am overshoot. One of many solutions to eliminate 
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the overshoot in step response is applying an auto-tuning stabilizer of PID structure, shown in Figure 6. PID 
controller is a generic control loop feedback method that will correct the error between a measured process 
variable and the desired input by calculating and give an output of correction that will alter the process 


accordingly. A PID controller has the general form: u(t) = K,e(t) + K; Í i e(t)dt + Ky = where K,,, K;, and 
Kq are proportional, integral, and derivative gains, respectively. The calculation involves three separate 
parameters; the proportional, the integral and derivative values. The weighted sum of these three actions is 


used to adjust the process for obtaining the desired output. The tuning rule based on process reaction curve 
methods or Ziegler-Nichols’s method [31], 
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Figure 6. Block diagram of the proposed auto-tuning stabilizer of PID structure 


2.4. Simulator of auto-tuning stabilizer of PID on transient generator under MATLAB/Simulink 

A linearized model of synchronous generator unit is programmed shown in Figure 7, using basic 
Simulink model components. As shown, the inputs are the excitation command voltage, Ep, stator winding 
voltages, Vae and Vje. The outputs are terminal voltage magnitude, V, terminal current, /, active and reactive 
powers at terminals, Fyen and Qgen, relative load angle, ô, slip, (w, — we), and accelerating torque, Tem. 

In Figure 8, a model of an excitation system is shown. As shown, the inputs are the excitation 
command voltage, V,ef, the supply voltage, Vsupp, obtained from output of auto-tuning stabilizer of PID 
structure, and terminal voltage of generator, V;. A model of an exciting system may be represented by the 
third order model. As shown, the voltage equation of the field winding and region of magnetization curve, S,, 
around the normal operation condition can be approximate by an exponential function, A,,exp?e*”*. 

A single-machine infinite bus system connected with exciter and auto-tuning stabilizer of PID 
structure can be easily constructed from the modules described in section 2, as shown in Figure 9. This can 
be accomplished by choosing specific number of generating unit and connecting them to an infinite bus. 
Once the auto-tuning stabilizer of PID structure is interfaced with the single-machine infinite bus system, an 
output of the network injects the exciter. And finally, the exciter will inject the transient generator. 

We can access GUI facilities of MATLAB to compose a software package of simulator for studying 
the effects of PSS and PID on transient synchronous generator. The synchronous generator and conventional 
PSS parameters are shown in Table 2. As an example of using GUIDE, the MATLAB GUI development 
environment, abilities, menu and plotting commands are employed in a script file or m-file to deliver 
interactive windows. These tools significantly reduce the process of designing and building GUIs. GUIDE 
automatically produces a MATLAB program file that controls how the GUI works [32]. 


Table 2. Generator, exciter, and PSS parameters 
Unit Parameters 
Generator V,-ated=18 KV Prateq=83 MW 1,=0.0048 pu xg=1.790 pu 
Xq=1.660 pu x/,=0.215 pu xg =0.355 pu x4=0.570 pu 
xq=0.275 pu x4 =0.275 pu Tg,=7.9 pu Tgo=0.410 pu 
Tjo=0.032 pu Tj4.=0.055 pu H=3.77 pu D,,=2 pu 
KA=50 pu V; = 1.0 + j0.0 S; = 0.8 + j0.6 
Infinite bus Tə + jXe=0.027 + j0.1 (high inductance) re + jx,=0.013 + j0.05 (low inductance) 
Exciter TA=0.06 Vrmax= | pu Vemin=-1 pu 
TE=0.052 KE=-0.0465 pu TF=1.0 pu 
KF=0.0832 AEx=0.0012 pu BEx=1.264 pu 
PSS K;=120 Ty=1.0 T1=0.024 
T2=0.024 T3=0.024 T4=0.002 
PID Controller K, =1.3 K;=0.9 Kg=0.2 
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To perform different operations and to visualize the different dynamic performances we have 
fulfilled and established a GUI under MATLAB. This Interface permits as to improve human-computer 
interaction, present control system from conventional PSS, fast view of the system regulation results and 
simulation, compute the system dynamic parameters, and exploit of the different PID controller, shown in 
Figure 10. 
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Figure 8. Simulink model diagram of exciting system 
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Figure 9. Simulink model diagram of generator connected with exciter and stabilizer of PID structure 


SIMULATOR OF PSS AND PID EFFECTS 
ON TRANSIENT SYNCHRONOUS GENERATOR 
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Figure 10. The window of inserting the inputs and displaying the variable assessments 


3. RESULTS AND ANALYSIS 

With the purpose of exploring the effect of auto-tuning stabilizer of PID Structure which comprise 
conventional PSS and PID controller, on transient synchronous generator, numeral simulations were carried 
out for several circumstances under different operating conditions, such as stable system and unstable system. 
The SMIB model is used and the values of delivered complex power S; and torque perturbation about 
operating point Tp are subject to change in pu. Figure 11 (a) and 11 (b) show Bode plots of GEP(s) with 
infinite buses of r, + jXe=0.027 + j0.1 and 7% + jx,=0.013 + j0.05, respectively; While Figure 12 (a) and 
12 (b) show Bode plots of PSS(s) with infinite buses of r, + jxe=0.027 + j0.1 and rm + jx,=0.013 + 40.05, 
respectively. The system operating points are at S; = 0.8 + j0.6 pu and V, = 1.0 pu, shown in Table 2. The 
gain of PSS(s) is DPSS=6. The phase plot of GEP(s) has two poles near 6.5 Hz or 40.8 rad/sec. From the 
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phase plot of PSS(s), we can see that the desired phase characteristic of the PSS ought to be lagging below 
40.8 rad/sec and leading above 40.8 rad/sec. There is no difference between two figures over desired 
stabilization frequency range centered around oscillation freq. of 40.8 rad/sec. The phase compensation can 
be designed using lag-lead and lead-lag units. The initial values of T} and T} are 1/40.8 or 0.024. And the 
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values of T, of lead-lag unit and T; of lag-lead unit might be 0.02 and 0.24, respectively. 
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Figure 11. Bode plots of GEP(s) in term of infinite bus’s value, (a) infinite bus: r, + Jx,=0.027 + j0.1, 
(b) infinite bus: re + jXe=0.013 + j0.05 
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Figure 12. Bode plots of PSS(s) in term of infinite bus’s value, (a) infinite bus: r, + jx,=0.027 + j0.1, 
(b) infinite bus: r, + jx,=0.013 + j0.05 


Figure 13, 14, and 15 show the plots of values of output voltage |V,|, delta angle 6, and active power 
Pg, respectively. As a result of a given value K, = 120 of wash-out network of PSS, the system response is 


stable on considerable perturbation, shown in those figures. When a LFO is appeared during unbalance 
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between power demand and available power at a period of time, for example, the instability becomes 
apparent in the system simulation shown in Figure 13 (b), 14 (b), and 15 (b), where a given value of wash-out 
of PSS is K, = 240 or heavy perturbation. Definitely, the shape of output voltage of unstable system is not a 
straight line but slightly wavy, shown in Figure 13 (b). Comparing both delta angles in Figure 14 (a) and 
14 (b), we can obtain significant difference; the shape of delta angle regarding the unstable system becomes 
out of under-damped response and has a big value of under-shoot, shown in Figure 14 (b). Figure 15 (a) 
shows a plot of active power’s variable that has a slightly slower of steady state response comparing to its 


counterpart in Figure 15 (b). 
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Figure 13. Bode plot of system responses in term of output voltage, (a) with considerable perturbation, 
(b) with heavy perturbation 
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Figure 14. Bode plot of system responses in term of delta angle, (a) with considerable perturbation, 
(b) with heavy perturbation 
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Figure 15. Bode plot of system responses in term of active power, (a) with considerable perturbation, (b) with 
heavy perturbation 


Figure 16 (a) and 16 (b) illustrate the results of delta angles and active powers regarding system 
response after using conventional PSS and auto-tuning stabilizer during instability condition, respectively. 
The performance of auto-tuning stabilizer of PID structure is better than conventional PSS in term of respon 
system during a heavy perturbation, shown in Figure 16 (b); i.e, that eliminates both over-shoot and under- 


shoot. 


Bulletin of Electr Eng & Inf, Vol. 10, No. 4, August 2021 : 2383 — 2395 


Bulletin of Electr Eng & Inf ISSN: 2302-9285 02393 


delta 
reoen 


— = 
SF SS _. 
= 








(b) 


Figure 16. System respons of delta angle and active power during a heavy perturbation; (a) with 
conventional PSS during a heavy perturbation, (b) with auto-tuning PSS of PID structure during a heavy 
perturbation 


4. CONCLUSION 

In this study, an auto-tuning conventional PSS of PID structure is additional controller to damped 
low frequency oscillation and to improve dynamic performance of the generating unit. The use of 
combination between conventional PSS and PID controller to study power system stability concepts 
regarding transient synchronous generator working at various operating conditions is investigated. Results 
from this study indicate that auto-tuning PSS of PID structure gives much better dynamic performance as 
compare to that of conventional PSS. To use a PID controller tuned by Ziegler-Nichols’s method will 
preserve the stability margin. Modeling of proposed auto-tuning stabilizer in Simulink environment provides 
an accurate result when compared to mathematical design approach. This experiment provided students of 
the electrical engineering study program with some realistic and challenging design experience and exposed 
them to a well-known power system stability problem regarding transient synchronous generator and its 
enhancement. The course survey indicated that the students were generally pleased with the design activities. 
Any one of the PSS design projects can be implemented for use with modifications of similar courses. A 
calculate algorithm can be added to obtain quantitative values regarding system response in term of a simple 
criteria. For a higher-level course, a PSS design for multiple-machine systems or the dual-input PSS design 
can be addressed. In addition, PSS designs for systems with multiple operating conditions can be counted. 
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